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The SUSY hybrid inflation model is found to dissipate radiation during the inflationary period. 
Analysis is made of parameter regimes in which these dissipative effects are significant. The scalar 
spectral index, its running, and the tensor-scalar ratio are computed in the entire parameter range 
of the model. A clear prediction for strong dissipative warm inflation is found for n$ — 1 — 0.98 
and a low tensor-scalar ratio much below 10 -6 . The strong dissipative warm inflation regime also 
is found to have no ^-problem and the field amplitude much below the Planck scale. As will be 
discussed, this has important theoretical implications in permitting a much wider variety of SUGRA 
extensions to the basic model. 
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I. INTRODUCTION 



The success of the inflationary paradigm has motivated in recent times more serious efforts in building realistic 
' particle physics models that incorporate cosmology 0,00. The objective of this sort of model building is to account 
for various cosmological features, central being inflation, but also leptogenesis, dark matter etc... to constrain the 
high energy properties of the model, and such that in the low energy regime the model reduces to the Standard 
Model. Supersymmetry (SUSY) plays a central role here, since aside from its attractive features for particle physics, 
it also allows stabilizing very flat scalar potentials, which are essential in inflation models due to density perturbation 
constraints. In this respect, a widely studied SUSY model of inflation has been the hybrid model, 

W = K 5($!$ 2 - /i 2 ) , (1) 

*-£h I where $i, $2 are a pair of charged fields 1 under some gauge group G, and S is the singlet which plays the role of the 
inflaton. 

O 1' An important feature about the model Eq. Q and its various embeddings into more realistic particle physics models 
00 is that the inflaton field generally interacts with other fields, with coupling strengths that can be fairly large. 
Even though the nonzero vacuum energy necessary to drive inflation will break SUSY, this underlying symmetry can 
stiilprotect the very flat inflaton potential from radiative corrections arising from these perturbatively large couplings 
0, 0- It has been observed in recent works that the effect of interactions of the inflaton with other fields 
does not simply affect the local contributions to the inflaton effective potential, but also induces temporally nonlocal 
■ terms in the inflaton evolution equation, that in the moderate to large perturbative regime yield sizable dissipative 
effects. Although SUSY cancels the large local quantum effects, for the dynamical problem the nonlocal quantum 
effects can not be canceled by SUSY. These dissipative effects in general can lead to warm inflationary regimes 0. 
Thus the conclusion of the works 0, 13 is that in general, models in which the inflaton has interactions with 
other fields with moderate to strong perturbative coupling, inflation divides into two different dynamical regimes, 
cold [Tol [Til IT^ ] and warm 0. This finding is very important, since these two types of inflationary dynamics are 
qualitatively much different. Thus one should expect different observable signatures in the two cases, as well as other 
theoretical differences in the treatment of inflation. 

The purpose of this paper is to apply these recent findings about dissipative dynamics during inflation to the SUSY 
hybrid model Eq. Q and to common extensions of this model. In particular, two models will be studied in this 
paper, Eq. and this model with a matter field A A coupled to it as 

W = /cS($i$ 2 - fJ 2 ) +.9$ 2 AA. (2) 



'Electronic address: mbg@ugr.es 
^Electronic address: ab@ph.ed.ac.uk 

1 One could also consider the fields $1,2 being gauge singlets. However, in this case the gravitino constraint on the reheating temperature 
translates into k being at most of the order of 10 — 5 H. That constraint is avoided when instead $1 2 are non-singlets. 
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The above is a toy model representing an example of how the basic hybrid model Eq. is embedded within a more 
complete particle physics model, in this case through the A fields. 

In this paper we will study inflation for both models Eqs. and We will show that in the above models both 
cold and warm inflation exist and we will determine the parameter regime for them. This will then explicitly verify 
the conclusions from the recent papers on dissipation 0,13, that showed both types of inflationary dynamics could 
exist. In both inflationary regimes, we will calculate the scalar spectral index ns — 1, and its running dns/dhik. 
With this information, we will then identify the qualitative and quantitative differences arising from the warm versus 
cold regimes. 

We emphasize here that the main objective of this paper is to determine in an explicit and well-known model how 
these two dynamically very different inflationary regimes emerge. This result departs radically from current wisdom, 
in which it is tacitly always assumed that all regimes in any model are cold inflationary. In light of our result, many 
other particle physics model building issues emerge in the warm inflationary regimes that we identify. Although we 
make some effort in this paper to discuss these issues, it is not the purpose here to dwell on them. In cold inflation, 
these issues have been the subject of many years of study, and likewise a complete understanding of similar issues for 
warm inflationary regimes will take focused effort in future work. 

The paper is organized as follows. Sect. [U] reviews the effective potential for the SUSY hybrid model, as well as 
some of the basic results based on this model for cold inflation. Sect. IIIII reviews the main results of dissipative 
dynamics during inflation based on 0, Q . Included here are the main formulas for density perturbations and the 
scalar spectral index in the presence of dissipation and a thermal bath. Subsects. IIII Al and IIII Bl compute the effect 
that dissipation has on inflaton evolution in the two models Eqs. and J2J respectively. The outcome of the analysis 
in these two subsections is a graph that divides the parameter space of the model into cold and warm inflationary 
regimes, and associated predictions for the scalar spectral index, its running and the tensor-scalar ratio. Sect. IIVI 
addresses some model building issues that emerge in the newly found warm inflationary regimes, such as the gravitino 
abundances and the constraint on the reheating temperature. Finally Sec. V states our conclusions. 



II. SUSY HYBRID MODEL 

We briefly review first some well known results about standard supersymmetric hybrid inflation [l3L IbH Il5| , in 
order to study later what are the main modifications introduced when taking into account the dissipative dynamics 
present during inflation 0, |^] . We consider the standard superpotential for the supersymmetric hybrid inflationary 
model, Eq. Q). Without taking into account SUSY breaking, the zero-energy global minimum in the model Eq. 

is located at the vacuum expectation values (VEV) S=0, <f>i = <f> 2 = M- If t ne g au g e group G of $i and $2 is 
identified with a GUT symmetry, the scale n would be the GUT symmetry breaking scale. On the other hand, for 
\S\ > [i, there is a local minimum at $1 = $2 = with potential energy given by the constant term k 2 /i . Inflation 
occurs while the system is located in this false vacuum. Here, the inflaton scalar field S and its fermionic partner- 
remain massless at tree level, while the scalars $1 and $2 combine into a pair of real scalar and pseudoscalar particles 
with mass = k 2 (|S'| 2 + /z 2 ), and another pair with mass m 2 _ = k 2 (|S'| 2 — /i 2 ); their fermionic superpartners are 
degenerate with mass mp = k\S\. 

Due to the splitting in the masses, there is a non- vanishing one loop radiative correction to the potential, AV, which 
provides the necessary slope and mass correction 2 to the inflaton potential to drive slow-roll inflation. In particular, 
the first and second derivatives of the effective inflaton potential are given respectively by: 



y^VVyV 3l (x 2 -l)(x 2 + l) , x 2 + l^ MW rM ... 

Av = ^e^r ln ? + xln l^i = ^^ Fl[xh (3) 



«VW f 2l (^ 2 -i)(x 2 + i) x 2 + i\ «VW r 



where M is the dimensionality 3 of the $1.2 representations, x = cf>s/(V2fj,), 4>S is the real part of the complex field S, 
and we are setting the imaginary components of all the fields to zero for simplicity. For large x, we have: F±[x] ~ 1/x, 
F 2 [x\~-1/ X 2 , 



2 We will not add any additional SUSY breaking mass term for the inflaton during inflation, so that its mass is given purely by the 
radiative corrections. Given that the fi scale is typically of the order of the GUT scale, SUSY breaking masses of the order of O(l) TeV 
are negligible unless the coupling k < O(10 — 5 ). 

3 We take J\f = 1 throughout this paper unless otherwise explicitly stated. 
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Standard (cold) "slow-roll" inflation is characterized by having small slow-roll parameter and rjji <C 1, with 



„„ = m ^ = c » - f ^y ^ . «i) 



ta/\2 2 / \2 2 / \ 2 » ro 

V \ K I urn p \ . r r o ft / Km p \ jV 



V \ 47T/Z 



where rap — 2.4 x 10 18 GeV is the reduced Planck mass. Therefore, during inflation the evolution of the inflaton field 
is well approximated by 

ZH4> S + V ~ . (7) 

At this point one should worry about supergravity (SUGRA) corrections to the inflaton potential. Generically, those 
give rise to mass corrections for the scalars during inflation of the order of 0(H 2 ), with tjh ~ 0(1), spoiling inflation 
(the so-called r\ problem [Hi. lltj]'). However, such corrections are not present if we take the minimal Kahler potential 
for the fields with a superpotential like Eqs.Q) and J2J, and this is the choice we adopt in this paper. Nevertheless, 
they will induce a quartic term (plus some higher order corrections) in the inflaton potential 0, Ea , with 

The quartic term dominates the inflationary dynamics when <f>s ~ nip, which happens for k ~ O(l). 

The values of the coupling k and the scale fj, consistent with the inflationary dynamics are obtained by demanding 
that (a) we have "enough" inflation (at least 60 e- folds), and (b) that the amplitude of the primordial spectrum 
generated by the inflaton vacuum fluctuations are in the range given by COBE observations. The former constraint 
gives the value of the inflaton field N e (~ 60) e-folds before the end of inflation, 



x N 



which is then used to evaluate the amplitude of the primordial curvature spectrum, 

3 



^ /2 "(£)(i)"Vi(^) GS) i^'v^fe) ■ (10) 

Therefore, using the COBE normalization [HHj P^ 2 = 5x 10" 5 at iV e ~ 60, we have 4 that /x ~ 2x 10" 3 m P ~ 5x 10 15 
GeV. 

Given that implicitly we are working with a SUGRA model, at most the VEV of the inflaton field could be of the 
order of the Planck scale, preferably below that scale. From Eq. @ one can see that <p> rnp for values of K > 0.8 and 
so are excluded 01- Moreover, taking into account the quartic SUGRA correction in Eq.©, the spectrum becomes 
blue tilted (rig > 1) pi [171 \18L W\ already for k ~ 0.05, which is not favored by the observational data on the spectral 
index from WMAP 5 , n s = 0.93 ± 0.03 " 



III. DISSIPATIVE DYNAMICS DURING INFLATION 



Dissipative effects can be important already during inflation, modifying the inflationary dynamics described by Eq. 
10. These are related to the quantum corrections in the effective potential of the background field. When neither 
the decay of the inflaton nor that of the fields coupled to the inflaton are kinematically allowed, loop corrections 
to the propagators are real, and they are absorbed into the renormalized masses and couplings, order by order in 
perturbation theory. On the other hand, when the inflaton (or the fields coupled to the inflaton) can decay into other 



4 Corrections to this estimation appear for small values of k < 0.01, for which ~ O(l). 

5 From the WMAP data only, the best fit value without running of the spectral index is ng = 0.99 ± 0.04. A recent analysis of 
the combined Lyman-a forest spectra + WMAP |2 J also gives nc = 0.99 ± 0.03 with no running, and ns = 0.959 ± 0.036 with 
dns/dlnk = -0.033 ± 0.025. However, a similar analysis in Ref. |2jj gives n s = 0.98 ± 0.02 with dn s /d\nk = -0.003 ± 0.01. 
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particles, the propagators in the loop have the standard Breit-Wigner form, with an imaginary contribution related 
to the decay rate T. Therefore, when computing the 1-loop effective potential for the inflaton field, the contributions 



associated to the decay rate lead to dissipative effects ja, |a| • In general this reflects itself in the form of temporally 
nonlocal terms in the inflaton evolution equation. Under certain approximations this translate into a simple effective 
friction term T5 in the equation of motion for the background inflaton field |y, |8j , 

j>s + {3H + r s )<p s + AV' = 0. (ff) 

The emergence of this friction term due to these underlying decay channels implies the dynamics of the system is 
such that part of the inflaton energy is dissipated into the lighter particles produced in the decays, i.e., into radiation 
p R , with 

p R + 4Hp R = ? s $s- ( 12 ) 

Although the basic idea of interactions leading to dissipative effects during inflation is generally valid, the above set 
of equations has strictly been derived in 6, 8] only in the adiabatic-Markovian limit, i.e., when the fields involved are 
moving slowly, which requires 

^<H<r, (13) 

<PS 

with r being the decay rate. The second inequality, H < V is also the condition for the radiation (decay products) to 
thermalize. 

Thus in general any inflation model could have two very distinct types of inflationary dynamics, which have 
been termed cold and warm 0, . The cold inflationary regime is synonymous with the standard inflation picture 
llOUlfl lT ^ . in which dissipative effects are completely ignored during the inflation period. On the other hand, in the 
warm inflationary regime dissipative effects play a significant role in the dynamics of the system. A rough quantitative 
measure that divides these two regimes is p]( 4 ps H, where p]( 4 > H is the warm inflation regime and p]( 4 ~ H 
is the cold inflation regime. This criteria is independent of thermalization, but if such were to occur, one sees this 
criteria basically amounts to the warm inflation regime corresponding to when T > H. This is easy to understand 
since the typical inflaton mass during inflation is ps H and so when T > H, thermal fluctuations of the inflaton 
field will become important. This criteria for entering the warm inflation regime turns out to require the dissipation 
of a very tiny fraction of the inflaton vacuum energy during inflation. For example, for inflation with vacuum (i.e. 
potential) energy at the GUT scale ~ I0 15_16 GeV, in order to produce radiation at the scale of the Hubble parameter, 
which is ps f0 10_11 GeV, it just requires dissipating one part in 10 20 of this vacuum energy density into radiation. 
Thus energetically not a very significant amount of radiation production is required to move into the warm inflation 
regime. In fact the levels are so small, and their eventual effects on density perturbations and inflaton evolution are 
so significant, that care must be taken to account for these effects in the analysis of any inflation models. 

The conditions for slow-roll inflation (cpg <C U, </>s <C Hips) are modified in the presence of the extra friction term 
Tg, and we have now: 

£ T = 77^. ( 14 ) 



(i + ry 

where r = T g/(3H), and eh , r/n are the slow-roll parameters without dissipation given in Eqs. (JjjJ and JSJ. In addition, 
when the friction term T5 depends on the value of the inflaton field, we can define a third slow-roll parameter 

(I + r) A 

with 

* = 3H*xt" (17) 

Similarly to the slow-roll regime without dissipation, when r/r < 1 , ex < 1 , and enr < 1, Eqs. (|f 1 [> and (| 1 2fl are well 
approximated by: 

; AV 1 

(18) 



3H 1 + r ' 



r 



PR ~ —^^--—.enV, (19) 
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and the number of e-folds is given by: 



AV 



7 (l + r)#. 



(20) 



Obviously, when T5 <C 1 we recover the standard "cold" hybrid inflation (CHI) scenario. 

The effect of the dissipative term is twofold: on one hand, dissipation of the vacuum energy into radiation acts 
as an extra friction term and slows down the motion of the inflaton field, so that inflation last longer. That means 
that when T 5 is non negligible, we would require in general smaller initial values of the inflaton field in order to have 
"enough" (at least 60 e-folds) inflation, Eq. (|20|l . On the other hand, fluctuations in the radiation background affect 
those of the inflaton field through the interactions, and this in turn will affect the primordial spectrum generated 
during inflation. Approximately, one can say that when T > H the fluctuations of the inflaton field are induced by 
the thermal fluctuations, instead of being vacuum fluctuations, with a spectrum proportional to the temperature of 
the thermal bath. We notice that having T > H does not necessarily require T5 > 3H. Dissipation may not be 
strong enough to alter the dynamics of the background inflaton field, but it can be enough even in the weak regime to 
affect its fluctuations, and therefore the spectrum. Depending on the different regimes, the spectrum of the inflaton 
fluctuations P\^ is given for cold inflation [24| , weak dissipative warm inflation (2^, |2(| and strong dissipative warm 
inflation |2]J respectively by 



T <H : 



To < H < T 



T S >H 



pl/2, 




Jj F S<t> \T=0: 



(21) 
(22) 
(23) 



with the amplitude of the primordial spectrum of the curvature perturbation given by: 



P. 



1/2 



■R 



H 



P, 



1/2 



3H 2 



AV 



(l + r)P t 



1/2 



(24) 



Given the different "thermal" origin of spectrum, the spectral index also changes with respect to the cold inflationary 
scenario mm sip, even in the weak dissipative warm inflation regime when the evolution of the inflaton field is 
practically unchanged. Again, for the different regimes, it is obtained: 



T<H : n s - 1 = -6e H + 2n H , 

17 3 1 

T s < H <T : 715-1 = — —e H + -rj H - -/3r , 

T S >H: ns -l = {-\e H + ^ m -- A f3 r )j^^ 



(25) 
(26) 

(27) 



In the latter case there could be appreciable departures from scale invariance in the spectrum; we notice again than 
in the strong dissipative case, slow-roll only demands e# <C (1 + r) 2 and rjn <C (1 + r) 2 , whereas the spectral index 
depends on the ratios £///(l +r) and T)h/(1 +r), which not necessarily are much smaller than 1. 

The question then is not whether there is dissipation during inflation, but whether this will affect the inflationary 
predictions. First, how large can T5 be in a realistic set-up. In the calculations in 0,0)0 a robust mechanism for 
dissipation during inflation has been identified. The basic interaction structure for this mechanism is 



r 1 2 ±2 2 



(28) 



where <j) is the inflaton field, x an d are additional fields to which the inflaton couples, and tpd are light fermions 
into which the scalar x-particles can decay m x > 2m^ d . This interaction structure can be identified in both models 
we are studying in this paper Eqs. |QJ and J5J. In the next two subsections, the dissipative properties of these two 
models will be computed based on the results on 0, Q and then the effect of this dissipation to inflation will be 
studied. 
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A. Decay into massive fermions: 

First consider the model Eq. JQ), which has only the minimal matter content. In this model, the scalar with the 
largest mass m+ can decay into its fermionic superpartner $ + and a massless inflatino, with decay rate 6 : 

where again xjq is the value of the inflaton field (normalized by fj,) N e efolds before the end of inflation. This decay 
rate is always smaller than the rate of expansion during inflation: 

and strictly speaking the adiabatic-Markovian approximation would not apply. Nevertheless, in order to get some 
numbers, let us proceed and estimate the dissipative coefficient and the amount of radiation produced. The former is 
given by: 







V2(^/4)(T + /m + ) 4>% 



64^1 + (r+/m+)y vi + 1 m+ 

~ YVi^J (x 2 N + l) 5 / 2fl ~ ~~2 V47r/ (31) 
which is always suppressed with respect to the expansion rate during inflation, H ~ Kfi 2 / (^/3mp): 

Nevertheless, the amount of "radiation" produced, i.e., the energy density dissipated from the inflaton, could be larger 
than if 4 , 



H± ~ 23^^V"W!UJ {—) (t> n + 1)*/* Fi[xn] ■ (33) 

Given that the ratio pa/H 4 goes like the inverse of x%, the ratio increases as the value of the inflaton field decreases 
during inflation. However, it only becomes larger than one for k ~ 0.1 toward the end of inflation, well after the 60 
c-folds before the end. We can conclude then that dissipation in this example is too weak to affect either the spectrum 
of the primordial perturbations or the dynamics of the inflaton field. 

B. Decay into massless fermions: 

As we have seen, dissipation through the decay into massive fermions does not have much effect during inflation. 
However, in more realistic models, one would expect the presence of other fields, which in principle are not directly 
relevant during inflation but can play a role during/after reheating. For example, fields coupled to either $1 or $2 
are massless during inflation, and become massive at the global minimum. Thus the model Eq. Pjl. where we have 
introduced a pair of matter fields, A and A, coupled to $2. Because they are massless during inflation, they do not 
contribute to the 1-loop effective potential, and radiative corrections are the same as in the previous case, with the 
slope and curvature of the effective potential given by Eqs. J2J and Q. But now the heaviest field with mass m + 
can decay into the massless fermionic partners of A and A, with decay rate: 

T+ = ^~ra + . (34) 

167T 



The interaction Lagrangian is given by C = — (k/v / 2)0+S<I > ^ 
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Since now there is no phase space suppression factor in the decay rate, we have T + oc <fis, and this can be much larger 
than the Hubble rate during inflation: 

£=V5£feW + l)V>. (35) 



H 16tt V p 

Having T + /H > 1, all the way up to the end of inflation, only requires g > 0.16 for k < 0.001 ( g > 0.01 for k = 0.5). 
This allows us to apply the adiabatic-Markovian limit in the effective equation of motion for the inflaton background 
field, Eq. <| 1 1 ft - with the dissipative coefficient given now by: 



S 2 UtJ \16tt 
and the ratio to the Hubble rate is given by: 



J N 



(1+^)V2 



p , (36) 



g 2 \ x" N mp 



'iH USy/Z-K \16ttJ (l + x 2 ,) 1 / 2 fjL 



(37) 



which behaves like Tg/(3H) oc x^ oc (fig, and so decreases during inflation. That is, the evolution of the inflaton 
field may change from being dominated by the friction term Tg to be dominated by the Hubble rate H . Whether the 
transition between these two regimes happens before or after 60 e-folds will depend on the value of the parameters of 
the model like k and g. The amount of "radiation" obtained through the dissipative term, is given by: 

PR „ 9 r , w p f o Xl 
H 4 ~ 2(l + r) 2 «V ' 1 ' 

which even when Tg < H could give rise to a thermal bath with T > H . In particular, we can have: 
(a) T s > 3H, and T > H {fa - -V^/T s ): 



(b) T s < H (fa - -Vt/QH)): 



PR 36V5 1 / mp\ 1 



PR 9 ( k \ Y 9 2 \ (rnp\ 7 J_ 

H 4 ~ 256V3tt V4tt/ U 67r / V M / x N 



(39) 



(40) 



The values of the couplings k and g for which we could have cold or warm inflation, and strong or weak dissipative 
dynamics, are plotted in Fig. QJ. In order to get the different regions, we have proceeded as follow: for each pair of 
values in the plane k — g, the value of the inflaton field at the end of inflation is determined. This is done in the cold 
and weak dissipative regimes by the condition 7 rjx = 1, Eci. (|15f) . In the strong dissipative regime inflation can end 
either with rjt = 1 or it may also happen that most of the vacuum energy is already transferred into radiation during 
inflation, and then inflation will end when pr ~ K 2 p 4 instead. In this case, whichever occurs first fixes the value of 
the inflaton field at the end of inflation. The value of the inflaton field at 60 e-folds of inflation is then obtained from 
Eq. H20|) . This in turn fixes the value of the dissipative coefficient Ts, Eq. I|36() . the temperature of thermal bath, 
Eq. (|38l) . and therefore the amplitude of the spectrum Eq. (|24H . The COBE normalization is then used to fix the 
value of the scale p. In order to match the expressions for the spectrum across the different regimes, we have used a 
simple expression with : 



pl/2 



iH 2 



AV 




(41) 



We can see in Fig. (QJ that the strong dissipative regime Ts > 3H requires large values of the couplings, n ~ g ~ 
O(l); for values k ~ g ~ 0.1 we are in the weak dissipative regime; and for values k ~ g ~ 0.01 we recover the cold 
inflationary scenario. Typically, for a fixed value of the scale p the amplitude of the spectrum in the strong dissipative 



7 The value of r/x becomes larger than 1 before the other two slow-roll parameters. 
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FIG. 1: Regions of cold (T < H), and warm (T > H) SUSY hybrid inflation in the ft — g plane. The warm inflation region is 
divided into the weak dissipative regime with T > H and Ts < 3H (lighter shaded region), and the strong dissipative regime 
with Ts > 3H (darker shaded region). Included are also the contour plots of constant /i, and the adiabatic-Markovian limit 
r+ = H. The black region on the right of the plot is excluded because <j>s > mp. In addition, when SUGRA corrections are 
taking into account, values to the left and down the wide dot-dashed line are excluded. 



regime would be larger than the one generated at zero T. The COBE normalization implies then a smaller value of 
the inflationary scale /i. For example, for k — g = 1 we have \x ~ 10 13 GeV, whereas pushing the coupling toward 
its perturbative limit, k = g — \ZAtt we would get /i~2x 10 10 GeV. On the other hand, going from the cold to the 
weak dissipative regime, the value of n only varies by a factor of 2 or 3, and it is still in the range of the GUT scale 
O(10 15 ) GeV. 

As mentioned before, the quartic term in the infiaton potential induced by SUGRA corrections becomes non 
negligible for not very small values of k. In the CHI scenario, the value of the infiaton field becomes larger than 
mp for k > 0.15, and consequently that region is excluded. The same constraint applies also in the weak dissipative 
scenario. However, in the strong dissipative regime, with Ts > 3H , the extra friction term keeps the value of the 
infiaton field below the Planck scale, and the constraint on k can be avoided. 

In Fig. we have compared the prediction for the spectral index of the scalar spectrum of perturbations in both 
the CHI scenario, and warm hybrid inflation (WHI). From the warm inflation scenario we can always recover the CHI 
prediction by taking g <C 1 ■ In standard SUSY GUT hybrid inflation, for small values of the coupling k the spectrum 
is practically scale invariant, it reaches a minimum around K ~ 0.01, and then rises due to SUGRA corrections up 
to positive values, which are disfavoured by WMAP results. But in the weak and the strong dissipative regime, due 
to the different origin of the spectrum, we get that the spectral index is still below 1 even for values of the coupling 
K > 0.01. This is especially true in the strong dissipative regime, where the dynamic is such that the infiaton field is 
well below the Planck scale and SUGRA corrections are negligible. In that regime the departure from scale invariance 
is within the observational value, with n$ — 1 — —0.022. 

In Fig. @ we have also compared the running of the spectral index, dns/dhik in both scenarios. The running 
although negative is much smaller than the value preferred by the WMAP data, dns/dlnk ~ — 0.0311q'q]J. Again, 
in the strong dissipative regime we can have larger values of the couplings k and g compatible with observations, but 
the predicted running of the spectral index is still small, with dns/d\nk ~ —4 x 10~ 4 . Nevertheless, it is not yet 
clear the statistical relevance of this result, which in any case would be finally confirmed or excluded by the Planck 
Satellite experiment. 

Concerning the primordial spectrum of tensor perturbations, as they do not couple strongly to the thermal back- 
ground, the amplitude is the same than in CHI, with 

p -J-fE-Y (AO) 

tensor ~ m p \2ir) ■ K ' 
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FIG. 2: Spectral index for cold SUSY hybrid inflation (solid line, CHI), and warm inflation (gray region, WHI). The weak 
dissipative regime (T > H but Ts < 2>H) is given by the darker gray region (triangle down);the strong dissipative regime 
(Ts > 3H) is given by the light gray region (triangle up). 
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FIG. 3: Running of the spectral index for cold SUSY hybrid inflation (solid line, CHI), and warm hybrid inflation (gray region, 
WHI). The weak dissipative regime (T > H but Ts < 3H) is given by the darker gray region (triangle down); the strong 
dissipative regime (Ts > 3H) is given by the light gray region (triangle up). 

As we have seen, strong dissipation (k ~ 0(1)) translates into a smaller inflationary scale /x compared to that of 
standard CHI, therefore a lower value of H and a smaller contribution of the gravitational waves to the spectrum. 
However, the same level of primordial tensor perturbations can be obtained decreasing the value of k and no dissipation. 
This can be seen in Fig. |0J, where we have plotted the prediction for the tensor-to-scalar ratio, defined as 



Rg - 



Pn 



(43) 



versus the prediction for the scalar spectral index ng — 1, for the different regimes (cold, weak dissipation, and strong 
dissipation). The smallest values of k correspond to a practically scale invariant spectrum. At present, from the 
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FIG. 4: Tensor-to-scalar ratio versus the predicted scalar spectral index ns — 1. The cold regime, T < H , is given by the dark 
grey areas; the weak dissipative regime, T > H but Ts < 3H, is given by the striped area; and the strong dissipative regime, 
Ts > 3H , is given by the light gray region. 



Cosmic Microwave Background polarization measurements the tensor-to-scalar ratio is poorly constraint (R g < 0.4), 
although in the future ratios as low as 10~ 6 could be probed (3^1 • Still, such a gravitational background can be 
achieved (but not larger) in this kind of models only in the cold or weak dissipative regime, and typically for values of 
the coupling k close to the maximum allowed by SUGRA corrections (blue tilted spectrum). On the other hand, in 
the strong dissipative regime we obtain a clear prediction that distinguishes this regimes from the others: no expected 
tensor signal, with R g < 10 -9 , and a red tilted spectrum with ns — 0.98. 



IV. REHEATING TEMPERATURE 



One of the main constraints on model building of inflationary supersymmetric particle physics models comes from 
the gravitino constraint on the reheating T after inflation. Gravitinos with a typical mass of the order of 0(1 — 100) 
TeV can be thermally produced during the radiation dominated era that follows inflation. If T is too high, we will have 
too many of them, and their subsequent decay will interfere with the predictions of the abundances of light elements 
at the time of big bang nucleosyntheses (BBN) j3^| . This puts an upper bound on the reheating Trh typically of the 

order of 10 9 GeV [33|, and a more recent analysis on BBN has lowered this bound to Trh ~ 10 7 GeV for a gravitino 
mass m 3 / 2 — 0(1) TeV IH^I ■ I n c °ld inflation, the reheating T after inflation can be well approximated by: 



T RH ~ — y/f^, (44) 
V n 9* ) 

where is the decay rate of the inflaton field, and <7* the effective number of degrees of freedom (typically of the 
order of O(100) in a SUSY model). Therefore, the reheating constraint translates in the inflaton not decaying too 
fast after inflation, which may implies some further constraint on the couplings 8 . 

In this letter we have just minimally extended the inflationary sector by adding a pair of extra matter fields A, 
A, with generic coupling g. This allows the decay of the <&i fields already during inflation, and the possibility of 
having a warm regime of hybrid inflation, when the couplings are not very small, say n,g > 0.1. The question then 
is whether this large coupling g could give rise to a large decay rate, and therefore a too large Trjj. From the S, $i 
sector we would have scalars and pseudoscalars with masses mj = \/2nu (plus a massless state if the minimum is 



This can be avoided if the gravitino abundance is diluted by the entropy produced during the late-decay of some other particle l3(i . 
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along the D-flat direction $1 = $2)- On the other hand, the A, A sector gives both scalars and fermions with a mass 
toa = \/2<?M- In order to avoid the decay of the inflaton into these fermions, it is enough to require re < 2g, which is 
just a mild constraint on the values of the couplings. 

In order to complete the transfer of the energy density into radiation after inflation, the model has to include the 
decay of the inflaton field into other lighter fields, with coupling hs and T$ — ms/ig/ (871"). In the strong dissipative 
regime (re > 0.1) we have fi ~ O(10 12 — 10 14 ) GeV, and demanding T < 10 9 GeV gives /ij<5x 10~ 6 . A similar 
constraint is obtained in the CHI scenario and the weak dissipative regime, where now ~ O(10 15 ) GeV, but 777,5 
decreases with re < 0.1. 

As a well-motivated example, which combines inflation with leptogenesis and light neutrino masses [3?1 13^|. the 
inflaton can decay into right handed (s)neutrino fields vm (i= family index). The decay proceed through the non- 
renormalizable coupling ^i^iVRiVm, with decay rate, 



J_ I M, 

8tt 



2 

m s , (45) 



where Mi is the RH (s)neutrino mass. In the CHI scenario, with u ~ O(10 15 ) GeV, re ~ 10~ 2 , and 777,5 — 10 13 
GeV, the gravitino constraint Trh < 10 9 GeV translates roughly into M; ~ 10~ 3 /j ~ O(10 12 )GeV . Those values 
are also consistent with baryogenesis and light neutrino masses 38]. This kind of scenario is also viable in the warm 
inflationary regime. Being consistent with the observed baryon asymmetry and the atmospheric neutrino oscillations 
does not directly constraint the value of re but the value of 7715 ~ 10 13 GeV. In the warm inflationary regime the value 
of the scale jj, required for successful inflation reduces as we moved into the strong dissipative regime, ms is of the 
order of 10 13 GeV for re ~ O(l), and the gravitino constraint gives now Mj ~ 10 _3 /j, ~ O(10 10 ) GeV. Therefore, a 
model of warm inflation and leptogenesis without the need of small couplings would be viable and compatible with 
observations, in the strong dissipative regime. 

Nevertheless, in the presence of a thermal bath already during inflation, one could worry about the value of T at the 
end of inflation, specially in the strong dissipative regime. It would be premature to impose the gravitino constraint 
directly on that temperature. Taking into account the decay of the inflaton field, the entropy production during the 
reheating phase can dilute the abundance of the gravitinos thermally produced at the end of inflaton [3^ . Roughly 
speaking, the entropy dilution factor would be 7 = Srh I 'S en< i ~ T en d/TjiH, where the subindex "RH" and "end" 
refers to the end of reheating and the end of inflation respectively. One should study in more detail the reheating 
phase after "warm" hybrid inflation before drawing any conclusion, taking into account in addition that production of 
gravitinos during reheating does not take place in a pure radiation dominated universe. Inflation will end before the 
vacuum energy has been completely dissipated into radiation, and the singlets may still oscillate around the global 
minimum, with their energy density on average behaving like matter. The initial production of gravitinos would 
proceed initially in a mixture of radiation and matter, but would be later diluted by the entropy produced by the 
decay of the singlets. 



V. CONCLUSION 



The key result of this paper is that the SUSY hybrid model, in particular Eq. J2|, has regimes of warm inflation. Up 
to now, it has been assumed that this model in all parameter regimes has only cold inflationary dynamics. However, 
Fig. Q firmly dispels this belief, as it shows that the parameter regime divides into regions of both warm and cold 
inflation. In light of this finding, the scalar spectral index, its running, and the tensor-scalar ratio have been computed 
in the entire parameter space of these models. We find a clean prediction for strong dissipative warm inflation with 
775 — 1 = 0.98 and a tensor-scalar ratio that is effectively zero. As shown in Fig. (0J, this prediction is very clearly 
separated from the cold results, which up to now have been the expected predictions from these models. Also these 
predictions for strong dissipative warm inflation are clearly separated from those of weak dissipative warm inflation. 

Theoretically the effects of dissipation in these models also present distinctive features. In particular, in the strong 
dissipative regime there is no 77-problem. Moreover, even for large coupling re ~ 1, the inflaton field amplitude is 
well below the Planck scale. A consequence of these features is that SUGRA corrections are insignificant. One 
interpretation of this would be that a much richer variety of supergravity extensions of the basic model Eqs. Q and 
© are permissible in comparison to the cold inflation case. This could have important model building applications, 
especially when identifying viable inflation models in low-energy limits of string theory. 

The main purpose of this paper was to highlight the dissipative dynamics inherent in this very popular SUSY 
hybrid model and then to outline the variety of new features this implies. There remains a great deal about our 
results that must be studied in further detail in future work. For example, as shown in Q the process of radiation 
production will also induce small temperature dependent corrections to the inflaton effective potential. These effects 
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will alter predictions for density perturbations. In fact, as shown in |29j temperature dependent effects could introduce 
qualitatively new features to the scalar power spectrum, such as oscillations. Thus a more accurate treatment of 
density perturbations and a thorough examination of their evolution is important to consider in future work. Along 
similar lines, a deeper issue is that of thermalization. In the analysis in this paper, we followed the results in 0,0, 
which treats thermalization based on some simple criteria. As stated in those works, a proper dynamical treatment 
of thermalization is still needed, and the consequences of such work could make significant changes in regards the 
predictions for density perturbations in some parameter regimes. However this paper has outlined the basic result 
that there are vast parameter regimes in this model, in which there is particle production during inflation, and thus 
in these regimes the statistical state of the system is substantially altered from the ground state. 

For the parameters values of the strong dissipative regime, reheating may start with a fraction of the vacuum 
energy already converted into background radiation. The deca y p roducts of the inflaton acquire plasma (tempera- 
ture dependent) masses wich will affect the reheating process |40|, kinematically blocking the inflaton decay until 
the temperature falls below the inflaton mass. Processes involving different particle production (thermal or out-of- 
cquilibrium) mechanisms during reheating should then be reexamined, such as production of RH (s)neutrinos and 
leptogenesis. In any case, reheating is completed in the warm inflation scenario through new decay channels different 
from those active during inflation. Inflation per se does not force the couplings entering in Eq. J3J) to be small as we 
have seen, and neither does reheating and the Trh constraint. 

One general result that can be taken away from this paper is that warm inflation regimes can be expected in SUSY 
inflationary models. There are many other models aside from the one studied in this paper in which we expect to 
find warm inflation regimes. An interesting case, that is worth mentioning are SUSY models which lead to monomial 
inflaton potentials. One important example of such a model is the next-to-minimal supersymmetric standard model 
(NMSSM) H3 , which has the superpotential 

W = \<!>H u H d - k& + h t Q 3 U^H u + • • • , (46) 

where $ is a singlet superfield, H u , Hd are the Higgs doblets, Q3 the third generation left-handed quarks and U§ the 
corresponding right handed up quark, and ht the top Yukawa coupling. Identifying the singlet field <I> with the inflaton, 
leads to an inflaton potential ~ k 2 4> 4 . In standard inflation, such a possibility for the NMSSM has not been of great 
interest, since it is well known for such a chaotic inflaton potential that the amplitude of the inflaton is greater than 
the Planck scale, (<f>) > mp, thus leading to the problem of an infinite number of unsuppressed higher dimensional 
contributions entering the potential. However in warm inflation, it is known 8] that monomial potentials like this one 
yield observationally consistent warm inflation for field amplitudes below the Planck scale (4>) < mp. Thus, in warm 
inflation such potentials have no trouble with higher dimensional contributions, and so are completely consistent. As 
such, this fact implies that NMSSM is a model, with no further modifications, that can support inflation. This is one 
of the simplest and may even be the minimal model that is consistent with the Standard Model and yields inflation. 
In future work we plan to do a detailed analysis of warm inflation in the NMSSM. 
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